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A Perturbation Theory for Light Diffraction from 
a Bigrating with Multiple Surf ace-Polariton Excitation 



N. E. Glass 
Department of Physics 
Naval Postgraduate School 
Monterey, CA 93943 



Abstract 

A perturbation theory has recently been developed for 
treating the diffration of light, with surface-polariton resonant 
excitation, from a bigrating surface. The theory is an 
approximation, to first order in the grating height, of an earlier 
theory based on the Rayleigh hypothesis. The perturbation theory 
holds for arbitrary polarization and for arbitrary plane of 
incidence with respect to the grating directions. It is, however, 
limited to treating only the possibility of two simultaneous 
resonant evanescent waves, whereas four are possible at normal 
incidence on a square grating, and three are possible on a 
rectangular grating. The present work generalizes the earlier 
work to allow for a four-fold (or three-fold) resonance. This 
extension also allows one now to determine the complex dispersion 
relation for surface-polaritons , for wavevectors at the 
intersection of two Brillouin zone boundary lines (not possible in 
the previous formulation) . 
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I. Introduction 



The study of electromagnetic wave scattering from corrugated 

metal surfaces, with resonant excitation of surface-plasmon- 

polaritons, has many important technological applications. One of 

these is the creation of electromagnetic radiation from an 

electron beam passing above a grating surface -- the Smith-Purcell 
Cl) , 

effect' . In many applications, the enormous enhancements in the 

absorptance at the surface and in the field intensity in the 

surface region, due to the grating coupling of an incident wave 

into a surface-polariton (or surface-electromagnetic-wave, SEW) , 

can be further enhanced by use of a bigrating (a surface periodic 

in two different directions) in place of a classical grating 

(periodic in one direction) . It was pointed out in some of the 

( 2 - 3 ) 

earlier theoretical studies on bigratings' ' that the resonant 

absorption of unpolarized light (e.g., of sun-light in solar 

absorbers) will be greater on a bigrating than on a classical 

grating. When the plane of incidence is normal to the groove 

direction, only the p-polarized component of the incident wave can 

excite the surface-polariton (which is itself a p-polarized 

surface-wave) ; the s-component does not couple to the surface-wave 

at all. By contrast, the bigrating can couple both the p- and the 

s- components of the incident wave, regardless of the orientation 

of the plane of incidence. Moreover, as demonstrated 

( 4 \ 

experimentally by Inagaki et al. , the bigrating can couple the 
incident wave simultaneously into two surface-polaritons , leading 
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to an extra absorptance, for incidence orientations in which the 
classical grating would only excite a single polariton. These 
same processes indicate mechanisms by which the Smith-Purcell 
radiation can also be enhanced by using a bigrating in place of a 
classical grating. 

The first theory of the Smith-Purcell effect that was both 

mathematically exact and physically realistic was the work of P.M. 
( 5 ) 

van den Berg' . This theory used a Green's function approach, 

f 6 ) 

proposed by di Francia^ , in which the component evanescent 

waves, which make up the Fourier expansion of the field of the 

electron current, are seen as exciting the various diffracted 

orders of the classical grating. 

One shortcoming of the van den Berg theory is that it 

considers the metal on which the grating is ruled to be a perfect 

conductor. This approximation precludes the possibility of 

accounting for the effects of the surface-polaritons and their 

( 7 ) 

resonant excitation. Recently, Chuang and Kong' ' developed a 
theory which considers a sinusoidal grating on silver, described 
by a complex dielectric function, and uses an extended boundary 
condition method to calculate the radiated energy. 

Similar techniques are in principle applicable to the 
treatment of Smith-Purcell radiation from a bigrating, but the 
numerical difficulties can become prohibitive: the computing time 
necessary to calculate the diffracted field, arising from each 
Fourier component in the incident field, will be much greater for 
the two-dimensional case, the bigrating, than for the 
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one-dimensional case, the classical grating. Thus perturbative 
techniques are desired. 

A non-perturbative theory of light diffraction from a 

bigrating, with sur f ace-polariton resonances, was developed by 

/ 8 \ 

Glass, Maradudin, and Celli' . This theory uses the vectorial 
form of Kirchhoff's integral to eliminate the field in the metal 
and thereby to cut in half the size of the matrix equation to be 
solved. There are nevertheless very large matrix equations to be 
solved with this method: from 98x98 up to 243x243 in Ref. 8. 

Thus, a perturbation theory was recently developed by Glass'' , 
which applies the coupled-mode theory of Glass, Weber, and 
Mills ^ to the exact theory of Glass, Maradudin, and Celli^^. 
This perturbation theory reduces the computations to a 4x4 matrix 
equation, by retaining, to first-order, the two amplitudes (p-and 
s- polarized components) of the specularly diffracted beam and two 
evanescent-wave amplitudes (p-components only) corresponding to 
the possible simultaneous resonant excitation of two surface- 
polaritons. These two surf ace-polaritons may be propagating in 
non-collinear directions (unlike the case for the classical 
grating, in which two polaritons can be excited only for normal 
incidence, and then are propagating in exactly opposing 
directions). This theory as implemented in Ref. 9 has one 
shortcoming: it cannot treat the case of more than two 
surf ace-polaritons excited simultaneously. For example, it cannot 
treat the cases of normal incidence on a square-grating when there 
are four surface-polaritons excited simultaneously or of oblique 
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incidence on a rectangular grating when there are three 
simultaneous excitations. This shortcoming could be an important 
one if one wishes to use these perturbative approaches in 
calculating the Smith-Purcell radiation. Fortunately, it is not a 
difficult problem to correct. A straight-forward extension of the 
theory of Ref. 9 will produce a 6x6 matrix equation for four 
resonant-wave (p-) amplitudes and the two specular beam 
amplitudes. It is the purpose of this report to present the 
results of the extension of the perturbation theory of Ref. 9 to 
treat the cases in which there are up to four surface excitations. 

II. Theory 

The physical geometry is as follows. The nominal flat 

surface is in the plane x 3 = 0. The surface profile is described 

by x 3 =C (x N ) [where x || =x 1 x 1 +x 2 x 2 ]. The region x >f ) is 

vacuum; x 3 <f(X||) is dielectric, described by the complex 

dielectric function e (w)=e R (w)+ie I (w) . The surface profile 

function is periodic in two directons. Light of frequency w and 

wavevector ll is incident from the vacuum (with angle of 
• • • • — ^ 

incidence 0). The projection of k onto the surface is k |( , 

A 

which makes an angle <#> with x^ . 

The electric field in the vacuum is 
3(x) = ^ 1 (w,ic || ) exp [ iic (| • x |( - ia o (u,k n )x 3 ] + 

+ 2 2 S (w,^g) exp[ii^-x (| + ia o (u,Kg)x 3 ] . (1) 

In the above equation we have: 
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( 2 ) 



% “ J ll + 3 ' 

where G is a reciprocal lattice vector defined by the bigrating 



periodicity, which for a simple square-lattice grating is 



-» 2tt a a 

G = ( m i X l + m 2 X 2 ) 



for irij=0 , ±1 , ±2 , . . . 



("a" is the period); and where 

^ A A 

k (| = k q = - (cos<#> x 1 + sin<#> x 2 )sin0 
Also appearing in Eq.(l) is 



a o (u '¥ 



. 2.2 v 2 , 1/2 

(w /c - Kg) ' , 

. , v 2 2 . 2 . 1/2 

i(Kg - o /c ) 7 



for < w 2 /c 2 

2 2 2 
for k£ > <j /c 

Lj 



(3) 

(4) 

(5a) 

(5b) 



The amplitudes of the incident and scattered fields are: 

3 a (!"k„) » B H + < J 3 X *»> B 1 



2 1 (o,)c ) = [k..+x 



and 



K 3 



( 6 ) 



r ( o,%) = [Kg - x 3 ^Kgy ] A ||(^%) + ( x 3 x Kg)A l( o,%) (7) 

The p- and s- components for the incident wave are expressed in 
terms of B|| and , and for the scattered waves in terms of 

the Rayleigh coefficients Ajj and A^ 

The Rayleigh coefficients are the unknowns in the problem. 

To find them one could expand the field in the metal in a Rayleigh 
expansion and then apply the boundary conditions at the grating 
surface (this is the Rayleigh method) . The method used in Ref. 8 
reduces the number of equations in half by employing the vectorial 
equivalent of the Kirchhoff integral to eliminate the field in the 
metal. The effective boundary condition that results then yields 
two infinite sets of linear algebraic equations for the unknown 
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coefficients. These are Eqs. (2.18a) and (2.18b) in Ref. 8 and 
equations (8a) and (8b) in Ref 9. 

The equations for the Rayleigh coefficients contain the 
integral : 

I(“|G) = SS d 2 x„ exp[-iG*x n -iaC (x„) ] , (8) 

c cell 

where a c is the area of a unit cell on the bigrating surface, over 
which the integral is performed. 

The perturbation theory begins by expanding the integrand of 
"I" to first order in £, thus yielding: 

I(a|G) = 5g' 0 - iaf(G) (9) 

where 



f (£) = SS d 2 x (| C (Xu ) exp (-iG* x„ ) . (10) 
c cell 

Under the assumption that the dimensionless parameter af is small, 
we can substitute equation (9) into the equations for the A's, and 
(assuming that C(0)=0) find the first-order results: 



A„ (G) 



a 



= l 



1 C(S-G') [ c 5g* a ii^') ■ a ^ lA l ( ^’ )] 



. 

G 



a.-*-* 

GG 

c-*-* 

GG 



6 2,o _i/3 5 r< G > i 



% 



B n" a G,o B i ] 



and 
A ± (G) 



- ia 33 Y [a 33' A n < 3 ' > + b 33’ A i< 2 ')] 



G' 






8 



^) 



P <? 



t a 5,o B u +b 3,o B J 



(Ha) 



(lib) 



Here the prime on the summation indicates that G'^G, and we have 
employed the notation A^ (G) for A||(w,i^) (similarly for A^) and 
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have used the following definitions: 



P-g = a(w,Kg) + a 0 (w,k N ) 


(12a) 


a ^, = a(G), Kg) - a o («,Kg.) 


(12b) 


a GG' = X 3* (K G X K G' } 


(12c) 


A A 




b ^' = 


( 12d) 


K-> K-3. 




^ G G 

C GG ' D (Ig' a(o,K^)a o (w,K^. ) 


( 12e) 


* 

< 

< 




e G K G k ll a(w,K^)a o (w,k || ) 


( 12 f) 


r o 2 2l 1/2 

a(o,K£) = &(o)- 2 - 


( i2g) 


L C J 





As in Ref. 9, we assume that the amplitudes of the two 
components of the specular beam, A||(0) and A^(0), can dominate the 
diffracted orders. Unlike the case considered in Ref. 9, we now 
consider the possibility of having four , not two, resonant 
evanescent waves. We may have four wavevectors in the surface 
plane, corresponding to four different reciprocal lattice vectors, 
which satisfy the resonance condition: 



%* = ^11 + = 5 (u) , 



sp 



for r=l , . . . , 4 



(13) 



where 5 s p(<i>) is a surface-polariton wavevector corresponding to 

the frequency of the incident light. To zeroth-order in the 

roughness (the flat surface limit), the locus of points 5 (w), 

sp 

for a given frequency, defines a circle. Hence, to zeroth-order 

the four vectors 5-3 lie on a circle, at the four-fold resonance. 

If any one, or all four, of the 5^ are at resonance, then for 

GIT 

that, or those, wavevectors we have to zeroth-order: 
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0 



(14) 



a(<j,Kg ) + e (w)a o ((j,K£ ) = 
r r 

We are therefore going to assume that six amplitudes can be 
larger than all the rest, these being A^O), A^(0), A (| (G^ , 
A||(G 2 ), A|| (G 3 ) , and A| ( (G 4 ) ; and for each of these we can write an 
explicit equation, using Eq.(ll), with the six important terms 
separated out from the sum on the right hand side. For A ( |(l) this 
procedure gives us: 



a 



A„(D = 



trr- 1 ( l ni-p)c lp A ||( p) + C(l)[c 10 A n (0) - a 10 A i( 0)] 

P = 2 ,3,4 

+ C(D[e 1 B || - a 10 BJ + { ]> ’ C ( 1-j ) jA,, ( j ) - &1 . A ± ( j ) ] 

j 

- ^ C (1-P) a n A , (P) 

^ X P ± J J (15) 

p=2,3,4 



To simplify the notation we replace the variable or subscript G ^ 
everywhere by j alone (where j=l,..,4 are the resonant beams and 
j=0 the specular beam) . The prime on the summation indicates to 
sum on all j^0,l,2,3,4 . The nonresonant terms have now been 

segregated in the brackets on the right hand side. The equation 
for A |( (2) is obtained from Eq. (15) by simply interchanging the 1 
and 2; and similarly for A |( (3) and A u ( 4 ) • 

The specular amplitudes are: 



a 



00 



A i (0) = ia 00 I n- r > a or A ll (r) - ^- B 



r=l 



p o 1 



+ l a 



00 



{Y C (~j ) t a oj A|| ( j ) +b o j A 1 



(j) ] + l C(-r)b 0r A 1 (r)} (16) 



r=l 
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(0, - C(-r)c 0r A ||( r, -5^22_B 



r=l 



+ 1- 



i— O' 

c oo . 



f(-jnc oj A || (j)-a oj A 1 (j)] C(-r)a 0r A l( r)} . 



(17) 



r=l 



The amplitude of a given nonresonant wave can also be 

expressed with Eq. (11) , but now we drop all the other nonresonant 

terms from the sum on the right hand side: 

4 

A ± (j) = ia jj{ } C(j“P) a j p A n (P) + C (j) [aj 0 A n (°) + bj 0 A 1 (0) ] 

P=1 



A n(j) = 



+ nj)[ajo B M + b jo B i] } 

ia . . /-I 

{ 2 C(j-P) c jp A „(P) + C(j)[c jo A || (0)-a jo A ± (0)] 

J 3 p=l 



(18) 



+ C (j) [e j B || -a j0 B 1 ] } . (19) 

We now follow the development given in Ref. (9) . We 
substitute Eqs.(18) and (19) for the nonresonant amplitudes back 
into Eq.(15) for the the resonant amplitude A(l); and for the 



denominator c 11 we use 



K, 



C 11 1 + a ( 1) a (1) 



a(l)a o (l)+K‘ 

a(l)a Q (l) 



( 20 ) 



The 



where now stands for and a(j) for a(w,K;*.) 

3 0O oi 

denominator c^ vanishes to zeroth-order at the K^-resonance. 
We multiply both sides of the resulting equation for A., (1) by 



a 2 (l)a 2 (l) -K^, in the form 

a 2 (l)a 2 (l)-K^ = [a(l)a Q (l)-K 2 ] [a ( 1 ) a Q ( 1) +K 2 ] 



( 21 ) 



on the right hand side and as 
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2 cj 2 - cj 2 

a 2 (l)a 2 (l)-Kj = e ^ ( — ~) . ( 22 ) 

c c 

on the left hand side. Here we have 

2 , 1+e . 2 t .2 

w r = (— )c K r , ( 23 ) 

which is the frequency of a surface-polariton at wavevector K 



on the flat surface. Finally, we note that in the equation which 

we now have for A ( | ( 1 ) , we are left with an A ( |(l) on the right 

. . . ~2 
hand side in a term proportional to f , which arose from the 

nonresonant terms. We bring this term over to the left side, and 

use it to renormalize the flat-surface polariton frequency: 

4 



- 2 2 
w = (J 

r r 



I 1 ' 



2 

C a nr 
rr r 



( 1+6 ) (J 2 K 2 
' ' r 



[ l C(r-p)f (p-r)a pp a 



2 

rp 



P=1 

(P*r) 



+ 



Y f(r-j)f (j-r) aj . 

j 



^1^+a 

C jj 




for r=l,2,3, or 4. 

Here we employ the definition: 

= a (r)a Q (r) [a (r)a Q (r) - K r 2 ] 

Next we go back to Eqs. (18) and (19) for the nonresonant 
terms and substitute them into Eqs. (16) and (17) for A^(0) 



(24) 



(25) 



and 



A„ (0) . 

At this point we have six linear algebraic equations for 
An (0) , A jl ( 0), and A |( (r) [r=l,...,4] in terms of B |( and ; 

expressed in matrix form these equations are: 
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(n 2 -n 1 2 ) 


L 12 


L 13 


L 14 


Mi 


N 1 




A||(D 


L 21 


(Q 2 -n 2 2 ) 


L 23 


L 24 


m 2 


N 2 




A m (2) 


L 31 


L 3 2 


(Q 2 -n 3 2 ) 


L 3 4 


m 3 


N 3 




A|| ( 3 ) 


L 41 


L 4 2 


L 4 3 


(n 2 -n 4 2 ) 


M 4 


N 4 


• 


A|| (4) 


°1 


°2 


°3 


°4 


(1+P) 


Q 




A„ (0) 


R 1 


R 2 


R 3 


R 4 


S 


(1+T) 




A ± (0) 



( U l~ Mf) B || - 


N 1 B 1 


( U 2- M 2 ) B II - 


N 2 B 1 


( U 3-M 3 ) B || - 


N 3 B 1 


(U4-M 4 ) B n - 


N 4 B 1 


(W+V-P ) B,! - 


QB ± 


a 


00 

p7 b i 


(X-S)B n - (T+- 



In this equation Q = v/ (2irc/a) = a/X ; fi r = w / (2Trc/a) ; while the 
other symbols, namely the capital letters L through X, are defined 
in the appendix. 

The problem now is to solve, numerically, this set of six 
equations for the six unknown coefficients, the A's. Then these 
six coefficients can be substituted back into Eqs. (18) and (19) to 
find the nonresonant coefficients. With all the coefficients thus 
at hand, we can then determine the scattered field vector- 
amplitude, 2 S (i^), from Eq. (7) — for each diffracted order and 
for each evanescent wave. We can thus obtain the reflectance, or 
efficiency, associated with each of the grating's diffracted 
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orders, as well as the total field and field-enhancement near the 
grating surface. 

III. Conclusion 

The analytic formulation of a perturbation theory has been 
presented here, in complete analogy to that developed in Ref. 9 , 
in order to treat the special case, of normal incidence on the 
bigrating surface, which was not treated in Ref. 9 . Numerical 
implementation of the theory, in direct analogy to that in Ref. 9 
is now being carried out. This implementation will include a 
comparison of the perturbation theory to the exact ,i.e., 
non-perturbative , theory, in order to test the limits of validity 
of the former. When the conditions at normal incidence are such 
as to excite simultaneously four surf ace-polaritons , the 
absorptance is expected to be enhanced approximately four-fold 
over that found when a single surface-polariton is resonantly 
excited. 

Again, in analogy to the discussion in Ref. 9, the results 
presented here can be used to calculate the dispersion relation 
for surface-polaritons on bigratings. We take the matrix on the 
left hand side of Eq. (26) , cross out the fifth and sixth rows and 
columns (corresponding to the specular beams) , take the 
determinant of the resulting 4x4 matrix, and set it equal to zero. 
The primed summation symbols that appear in the definitions of the 
matrix elements must first be redefined to reinclude the specular 
beams, j=0. Solution of the determinantal equation by fixing 
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and solving for complex w will then yield both the dispersion 
and the life-time of the surface-polariton. The present 
formulation, unlike that in Ref. 9, will enable us to treat those 
points in the )c |( -plane where two Brillouin zone boundary lines 
intersect, that is, where there are four coupled surface 
polaritons. Work on numerically solving the dispersion relation 
is just getting started. 
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Appendix 

Eq. (26) employs the following definitions for the matrix 
elements . 



L = h if (r-p) c 
rp r L rp 



- Y C (r-j)C (j-p) (f 




(A-l) 



j 



4 




pp rp pO 



p=l 
(* r) 




(A-2) 



j 



14 



= -h 



if ( r ) a r o - l n r -P)nP)“pp a rp b pO 



P=1 

r) 



- y f (r-j)f (j) (f -a. n +a .b.^a.. 
L r -| -jo rj -jO' 



-a 



00 



'00 L 



iC(-r)c Qr + j C(p-r)f(-p)a a 



p=l 
(* r) 



pp pr Op 



- y C (j-r)f (-j) (f n .c. -a . a . ) a . . 
L VJ '** v J/ v Oj jr rj jO' 

j 



H7t[ I ?<P>f<-e>“pp a op + Y ? <3)C (-j) ( f 0 j c jo +a 0j )a jjj 

p = l j 



-^[1 

00 L 



f(P)C(-p)a pp a op b p0 + C(j)C(-j)(f oj -b jo )a jo . 



P=1 



= -a 



00 



if ( -r ) a 0r " l C(p-r)f (-p)a a 



P=1 
(* r) 



pp pr po 



- y C (j-r)f (-j) (f • a^.+a. b._,)a.. 
L * VJ ' * v J ' v jr Oj jr jO' 

j 



“oo[ l 



C(P)f(-p)a pp a p 0 b p0 + l f(j)C(-j) (f jo -b jo )a oj a j 



P-1 

4 



“oo[ I 



f(P)C(-p)a pp b p ‘ 



P-1 



+ Y C( 3 )f <-□) ( ^ )a 

-1 IID J JJJ 



= h r [iC(r)g r - J* £ (r-j)f(j) 



T . C . 

D rj 



a 



00 



'00 . 
D 



I f (j)f(-j)T j C oj 



(A- 3 ) 

(A-4) 

(A-5) 

jj. 

(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A-10) 



(A-ll) 



w 



(A-12) 



X 



a 



C 



oo e o 

oo^O 



“ooF C(i>n-j)T j a oj 

j 



(A-13) 



In the above equations the following additional definitions have 

been introduced: 
c . 



f . = 
mj 

f . = 




g . = c . - e . 

y D DO 3 



h = - 
r 



t a c 
r rr 



eoi 



t . 
D 



q .a . ./c. . 
y D DD DD 



1 



1 2 



27rc/a„ 



( A-14a) 

(A-14b) 
( A-15 ) 

( A— 16 ) 
(A-17) 



The terms a... b... c... e. were defined in the text, Eq.(12). 
They are dimensionless, as are f^j and g^ . The summations 2’ 
are, as in the text, over all j except j=0, 1, 2 , 3 , and 4 . 
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